Optimal Control of Quantum Rings by Terahertz Laser Pulses 
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Complete control of single-electron states in a two-dimensional semiconductor quantum-ring model 
is established, opening a path into coherent laser-driven single-gate qubits. The control scheme is 
developed in the framework of optimal control theory for laser pulses of two-component polarization. 
In terms of pulse lengths and target-state occupations, the scheme is shown to be superior to 
conventional control methods that exploit Rabi oscillations generated by uniform circularly polarized 
pulses. Current-carrying states in a quantum ring can be used to manipulate a two-level subsystem 
at the ring center. Combining our results, we propose a realistic approach to construct a laser-driven 
single-gate qubit that has switching times in the terahertz regime. 
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In recent years there has been wide interest in quan- 
tum control of nanoscale systems. One of the main moti- 
vations behind these studies arises from the possibiUties 
of using tailored laser-pulse sequences for logic opera- 
tions [l|. Semiconductor quantum dots and quantum 
rings (QRs) [5] are likely to play an important role in 
these far-reaching developments. Their atomlike prop- 
erties together with a high flexibility in size and shape 
construct an ideal playground for quantum control. 

A fundamental question in laser control of a general 
iV-level quantum system is controllability, i.e., if the con- 
trol target such as a certain eigenstate can be reached 
even in principle. For example, full population trans- 
fer in a parabolic quantum dot is not possible (in dipole 
approximation) [3]. The first steps towards laser con- 
trol of QRs were recently taken by Pershin and Picr- 
marocchi who theoretically analyzed the formation 
of currents through population transfer in narrow QRs 
subjected to a circularly polarized laser pulses. Similar 
approach has been used also to generate ring currents in 
circular biomolecules such as Mg~-porphyrin 5]. Quan- 
tum optimal control theory (OCT) [6] is a powerful tool 
to find optimal laser pulses for controlling a quantum sys- 
tem. The iterative scheme developed within OCT [7] con- 
verges monotonically to an optimal laser pulse for reach- 
ing the prescribed target, such as a desired final quantum 
state, at the end of the pulse. 

In this Letter we apply OCT to semiconductor QRs of 
finite ring width. We construct optimal two-component 
laser pulses that drive the QR from a given initial state to 
any predefined target state. These terahertz pulses gen- 
erate the desired transitions in significantly shorter times 
and higher accuracies than previously used circularly po- 
larized continuous waves (CWs) of finite lengths. Finally 
we sketch how the full control of current-carrying states 
in QRs enables the construction of a coherent laser-driven 
single- gate qubit. 

The time-evolution of our system is described by the 




FIG. 1: (a) Shape of the external confining potential for a 
quantum ring and an example of a circularly polarized laser 
field, (b) Energy-level spectrum of a quantum ring. The 
transitions are allowed along the dashed line so that Al — ±1. 



time-dependent Schrodinger equation 

ifi^*(r, t) = [Ho - fie{t)] t), (1) 

where e{t) = (exit), £y(t)) is the two-component laser 
field propagating in z direction. The interaction between 
the field and the electron is modeled in the dipole approx- 
imation (length gauge) with the dipole operator fi ^ —er. 
The static effective Hamiltonian for the semiconductor 
QR located on the xy plane is written as 
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where the potential part with r"^ — + y'^ consists of 
a parabolic confinement and a Gaussian peak located at 
the center. Together they determine the shape of a two- 
dimensional QR which is visualized in Fig. [TJ We choose 
the parameter values Swo = 10 meV, Vb = 200 meV, 
and c? = 10 nm in order to follow the energy and length 
scales of QR experiments [2] . The parameters yield a ring 
radius of ro k, 22 nm. The effective mass of the electron 
in GaAs semiconductor medium is to* = 0.067me. Using 
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the dielectric constant k = 12.7eo for GaAs, the effective 
atomic units [marked as (a.u.) in the results below] scale 
as Ha* — (m*/K^)Ha w 11.30 meV, a*^ — [m* / nja^ « 
10.03 nm, and = ft/Ha* w 58.23 fs. 

Figure ^ih) shows the energy-level spectrum for nine 
lowest states ■(/'j, where I is the angular momentum and 
index j — |^| + 1 enumerates the levels in energy. Due 
to the finite width of the ring, there are other radial 
bands which, however, appear at higher energies, e.g., 
the second band appears at E = 5.03. The selection 
rules allow transitions only between states with consec- 
utive angular momenta I = ...,—2,-1,0,1,2,..., i.e., 
along the dashed line marked in the figure. The transi- 
tion probabilities between the levels can be calculated 
as Pjk = I (V'fclAIV'i) For the first levels shown in 
Fig. [Tfb) we get relative transition probabilities P12 = 1, 
P23 = 1.04, P34 = 1.12, and P45 = 1.24. Hence, the 
two-dimensional nature of the ring is reflected also in the 
varying Pj^: in one-dimensional rings, the dipole matrix 
elements are constant along the transition lines [1]. 

We apply OCT to find optimal laser pulses for popula- 
tion transfer from the initial state = "^{t = Q) to the 
target state $/, so that both of these states are prese- 
lected from the QR eigenstates "(/ij . In OCT formalism we 
thus use the projection operator |$/)($/| as the target 
operator, whose expectation value is maximized at the 
end of the pulse (t = T). This corresponds to maximiz- 
ing the overlap 

Ji[*] = {-^{n^f) ($/|*(T)) = I (*(T)|<i>^) \\ (3) 

We also require that the fluence (time-integrated inten- 
sity) of the laser pulse is as small as possible by minimiz- 
ing the functional 
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dta(t) [el{t) + el(t)] 



(4) 



where the predefined function a(t) acts as a penalty fac- 
tor, which allows us to impose constraints on the enve- 
lope function of the laser pulse Q. We can apply, e.g., 
sinusoidal envelope functions as demonstrated in Ref. M- 
However, in the examples below we have set a = 1 (rect- 
angular pulses) in order to simplify the comparison be- 
tween OCT and CW approaches. As the second con- 
straint, the wave function and its complex conjugate need 
to satisfy the time-dependent Schrodinger equation. This 
condition is expressed in a functional form (in a.u.) as 



Me,^,x] = -2Im 



dt 



(5) 

where x(t) is the Lagrange multiplier. Now, the mini- 
mum of the total Lagrange functional J = Ji + J2 + J3 
can be determined by setting the variations with respect 
to e, and x independently to zero. This yields a set 




FIG. 2: (Color online) (a) Components €^(1) (blue) and ^y{t) 
(red) of the optimized (solid lines) and uniform laser pulses 
(dashed lines) for |l) —>■ |2) transition, (b) Spectrum of the 
optimized field. The dashed lines correspond to eigenfrequen- 
cies (jJi2 and 0^23. 



of control equations 

idt-i'it) = 77*(t), *(0) = (6) 
idtxit) = Exit), xiT) = $/ (<i>/|*(P)) , (7) 

e(t) ^ (xitM'fit)), (8) 

a 

which can be solved iteratively in order to find the opti- 
mized field e{t) at the end of the procedure |7|. For the 
forward and backward propagation of '^{t) and 
apply a rapidly converging scheme discussed in Refs.[3,B 
This computational scheme has been recently included in 
the octopus code [lo| . 

The advantage of OCT becomes obvious already in 
the simplest transition from ■01 to 'ip2 (denoted below 
as |l) — *■ |2)). Figure llja) shows the x and y compo- 
nents of the optimized pulse when the pulse length is 
fixed to T = 100 (~ 5.8 ps) and the penalty factor is 
a — 1. In comparison, the dashed lines show the CW: 
ecw(i) = A[cos{ujt)x ±sm{u)t)y], where ± denotes a± 
polarization, and A and to are the amplitude and fre- 
quency of the field, respectively. Here we choose uj = uji2 
and A — ^12 = 7r//ii2T, where VLr is the Rabi fre- 
quency. These parameters correspond to the 7r-pulse con- 
dition. The CW leads to an occupation of 0.99 in the sec- 
ond state which, as demonstrated below, is insufficient in 
terms of the controllability of the QR. In contrast, the op- 
timized pulse produces a value 0.9998 for the same pulse 
length. The underlying reason for the difference becomes 
evident from the spectrum of the optimized pulse shown 
in Fig. [2l[b). Namely, the pulse contains a considerable 
fraction of the next eigenfrequency W23, preventing the 
population flow into the third state. If the system is sub- 
jected to a CW instead, the occupation loss is inevitable, 
and becomes pronounced at shorter pulse lengths when 
the amplitudes are correspondingly higher. 

In Fig. [3] the yields of OCT and CW approaches for 
the |l) 1 2) transition are compared for different pulse 
lengths in a logarithmic scale. Considering the required 
pulse lengths for desired occupation accuracies, OCT is 
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FIG. 3: (Color online) Maximum occupation of the target 
state in transition |l) |2) as a function of the pulse length. 
The blue open circles correspond to continuous waves and red 
filled circles to the optimal-control result. The insets show 
the densities |^'(T = 100) ^ when the corresponding achieved 
occupations are 0.99 and 0.9998 for these pulse types, respec- 
tively. 

clearly superior to the CW approach. The relation be- 
tween the pulse lengths of these methods Tcw/T'oct for 
a certain occupation always doubles when the accuracy 
is increased by an order of magnitude. For example, the 
occupation of 0.99 requires Tew ^ 100 or Tqct ^ 25, 
whereas for an accuracy of 0.999 the CW needs to be 
increased up to Tew ^ 300 and the optimized pulse only 
to TocT ~ 35. The main benefit in very high occupa- 
tions achieved by short optimized pulses is the ability 
to perform a large number of successive operations (see 
below). In single transitions we can qualitatively assess 
the required occupation accuracies by checking the final 
densities of the propagated wave function. The insets of 
Fig. [3] reveal that occupations close to 0.999 are required 
for the target states in order to preserve a reasonable 
circular symmetry of the wave function. 

The benefits of OCT become even more pronounced 
in multilevel transitions. As shown above, there is only 
a single excitation path in the lowest radial band of our 
QR. The path is subjected to strict selection rules that 
allow transitions only between consecutive levels with 
= ±1. Hence, in the CW scheme the intermediate 
states require full population before occupation of the 
final state. Eventually this leads to very slow transi- 
tion processes sensitive to external perturbations of the 
real device. In OCT approach instead, full population of 
the intermediate states is not needed, since the compo- 
nents in the optimized pulse drive the occupation contin- 
uously toward the target state, still obeying the selection 
rules. For example, in |1^ — > |4^ transition the maximum 
occupation required for the intermediate states |2) and 
|3) is 60 . . . 80%, depending on the pulse length (here 
T = 50...200). 

The chirality of a QR state can be changed in multi- 
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FIG. 4: (Color online) Schematic picture of transitions from 
I = —1 to / = 1 (a) and from I = —2 to I = 2 (b) (upper 
panel), optimized fields for these transitions (middle panel), 
and the occupations of the states (lower panel). 

level transitions from I to —I (or vice versa). Figure [3] 
shows schematic pictures (upper panel) of transitions 
from to 1/^2 i^) ^iid from -0^^ to ■03 (b). The op- 
timized pulses given in the middle panel of Fig. [4] are ap- 
proximately symmetric around the midpoints at t = T/2, 
and at this point the polarization direction of the field 
changes from a+ to a-. Hence, a backward process can 
be obtained by applying a pulse where the x and y com- 
ponents are exchanged. 

The intermediate states, which in this case contain the 
ground state with / = 0, require only partial occupation 
in the transition process. The intermediate occupations 
shown in the lower panel of Fig. [4] reach a maximum of 

60% in both ^ = -1 ^ 1 and I = -2 ^ 2 transi- 
tions. Similar to the pulse shapes, the occupation curves 
are approximately symmetric around t = T/2. The fi- 
nal occupations of the target states are 0.998 and 0.996, 
respectively. 

In Fig. Owe plot the time-dependent electron densities 
|5'(r,t)p (colorscale) and the current densities j(r, t) = 
ih{^>y^* - ^*V^)/2m* (arrows) at different times in 
the transition from I — — 1 to ^ — 1, corresponding to 
Fig. [DJa). The circularly symmetric initial state in (a) 
carries a counterclockwise current of / « 0.11 //A. Dur- 
ing the propagation, the density and current localize first 
to a single maximum (b). At the midpoint (c) there are 
two density maxima along the ring and opposite currents 
between them along the two paths, arising from the par- 
tial occupations of ^ = ±1. The net current is zero at this 
point. As expected, the final state reached at T = 55 (d) 
has a symmetric clockwise current. 

The circular currents in QRs induce magnetic fields 
through the ring in accordance with the Biot-Savart law. 
The current-fiipping processes correspond to changes in 
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FIG. 5: (Color online) Time-dependent electron density 
|\t(r, i)p [light shade (yellow), low; black, high] and the cur- 
rent j(r, t) (arrows) at different times in the transition from 
I — —1 to I = 1 driven by the optimized pulse [see Fig. IHa)]. 



magnetic fields as Bz ^ —Bz at the center of the ring, 
where B = /zo//2ro « 3/xT for a single electron occu- 
pying the state I = ±1. If needed, higher fields can be 
obtained by occupying higher states, by decreasing the 
ring size, or by changing the semiconductor host mate- 
rial. 

Combining the results presented above, we suggest the 
following single-qubit gate construction. First, the accu- 
rate experimental shape of the QR and thus the static 
Hamiltonian Q have to be determined in order to al- 
low for a full control. This can be done, e.g., by tuning 
the potential parameters to reproduce the single-electron 
spectrum obtained in a transport measurement [ll[ . Sec- 
ondly, the QR is controlled by a sequential master pulse 
achieving first the initial excitation from I = to I = Im 
(qubit initialization) and then performing ^ — 
operations at desired times so that the opposite oper- 
ation is obtained by exchanging e^ and ey components 
of the current-flipping pulse. The magnetic field induced 
at the ring center changes the spin state (up/down) of 
an applicable subsystem, e.g., a magnetic particle 4] or 
an attached quantum dot. This is feasible within our 
picosecond-scale gate operation times, since the char- 
acteristic relaxation times are of the order of nanosec- 
onds :4|, and the decoherence times in quantum dots (of 
similar physical parameters) have been measured to be on 
the millisecond scale 



regime is routinely reached by, e.g., quantum cascade 
lasers, and the technology of pulse refinement is under 
rapid development [isj . 

In conclusion, we have shown that a complete con- 
trol of realistic single-electron quantum rings can be ob- 
tained using (quantum) optimal control theory. The op- 
timized pulses perform the desired operations at signifi- 
cantly shorter times and better accuracies than the con- 
tinuous waves, and in multilevel transitions they are not 
restricted by the need to fully populate the intermediate 
states. This enables a simple construction of a pulse that 
coherently flips the spin of a subsystem at the center of 
the ring at preselected times. 
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12|. Actual quantum computation 



would still require a two-qubit gate and a stable read- 
out scheme, but these developments, as well as detailed 
experimental considerations are beyond the scope of this 
study. However, we point out that terahertz frequency 
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